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ABSTRACT 



In this paper, we propose a complex approach to evaluate a function sum of two non- 
commuting non Hermitian operators. Then, it is proposed an explicit expansion of the 
evolution operator in the case of the neutral K meson system influenced by an external 
interaction. Then, the importance of the procedure is pointed out to consider the alge- 
braic expansion of the time evolution operator when ever the dynamics decouples the 
internal transitions and center of mass motion. 
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I. Introduction 



The temporal evolution of metastable systems is governed by a non Hermitian Hamil- 
tonian with non orthogonal eigenvectors corresponding to complex eigenvalues. The 
problem to determine an elegant and compact form for the evolution operator 

U{t) = exp[-iHt] where H = Ho + V (1) 

is connected with the more general issue to express explicitly an arbitrary function of 
the sum of two noncommuting matrix operators. One of the most tantalizing method 
to evaluate this matrix function involves the subtleties of complex analysis and it was 
already been developed in the particular case of Hermitian operators [ 1 ] . The purpose 
of this paper is to extend this method to non Hermitian operators. The starting point 
is the generalized Cauchy's formula [ 2 ] 

f(A) = ^ / f{z)Q{z)dz , (2) 
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here, the integral is extended over a contour 7 in the complex z plane which encloses all 
eigenvalues of A. It is then possible to obtain an integral expression of f(A) in terms 
of the well-known resolvent operator 

g A (z) = (zl-A)- 1 . (3) 

If A and B are two non Hermitian and noncommuting operators, the resolvent operators 
of A and of the sum A + B are given, respectively, by 

g A (z) = (zl-A)- 1 , g(z) = [zl-(A + B)]- 1 . (4) 

In the convergence region of the geometric series 



the following expansion for Q(z) holds 



B ' " 



Q ( Z ) zX _ A Z X-(A+B) ^T-^XJ _ 
zX-A n=0 v 7 
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which will be useful later. The right and left eigenvectors of the operator A are defined 
by the relations 

A\$i) =A^> , (^|A = A,(^| . (7) 

Contrary to the case where A is Hermitian, the sets and although com- 

plete, they are not orthogonal and ^ However, the following relation 

{9 i \$ j ) = (9 i \Q i )6 ij (8) 

holds, and therefore it is possible to generalize the completeness relation using the 
following decomposition of unity 

V <*«i**> ' { ) 

We propose to give a spectral expansion of the function f(A+B) in terms of their relative 
eigenvalues. If V is a closed contour enclosing the whole spectrum of the operator A + B, 
then we have 

f(A + B) = ^-J f{z)Q{z)dz (10) 



and from Eq. (6) 



1 °° r 

f{A + B) = 2ri E / [^aBTGa] dz . (11) 
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The matrix elements will be obtained by 

i r 00 
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iJ^|$ 2 ) l = (12) 
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^ 0{ tr } ( 1 w^T s 7*^r (*^^> s|$2) x 

^ J dzf(z) [(z - Ai)" 1 ^ - A.J" 1 ...^ - A^.J" 1 ^ - As)" 1 ] 



where we have used the relation (9) and we have defined {i, n} = {ii, i 2 , i n -i}- The 
indices ik run through the whole set of the eigenvectors as usual, whereas and ($2) 
are fixed. If we introduce the following function 



F{z) = f{z) [{z - X 1 )(z - K)...{z - X^z - A 2 )] 



-1 



(13) 



we can write 

{VilfiA + B)^) 

n=0 {i,n} 



^~^J B ^-^J- ,> m ^ 1 dzF{z) ' 



(14) 



which generalize the result of the previous paper [ 1 ] . Supposing that the eigenvalues 
of A are enclosed within T, i.e. all the singularities of the function F(z) are inside the 
integral contour, it is then possible to apply the theorem of the residues. If we denote 
TZ(Xi) as the residue of F(z) at the pole z = Ai, the matrix elements in the Eq. (14) 
can be rewritten as 



n=0{i,n} V ll{ 11 ' \ *2l »2/ 



In — 1 / \ — 1 



■B\$ 2 ) x 



x 



n-l 



^(A!) + ^^(A,J+^(A 2 ) 



(15) 



If some eigenvalues are degenerate, the general expression of F(z) is 

F(z) = f(z) [(z - Ai) mi (z - K )<»*i ...(z - K^)^ (z - A 2 )-] - 1 (16) 

where every exponent r»i is the degeneracy order of the respective eigenvalues Aj. Finally 
we want to stress that the set {m^} (i £ {1, n — 1}) depends on the particular {i, n} 
selected. It is worth noting that these results recover the formulae already present in 
literature in the limiting case of the Hermitian matrices, and they result generally more 
straightforward than the usual algebraic methods [3]. The use of these results can be 
displayed in the practical example of the evolution operator. 



II. The evolution operator of the neutral kaon system 

The previous results let us make a decisive step toward a complete understanding of the 
controversial results about the dynamical behaviour of a decaying system described by 
the vector state \^(t)). Its time evolution can be written by means of an operator U: 

which can be expressed in the well-known exponential form, (using units h = 1) 

U(t) = exp[-iHt] . (18) 

Although the Hamiltonian of a sensible quantum system is expected to be a Hermitian 
operator, under suitable conditions we may recover the time evolution according to an 
effective non Hermitian Hamiltonian like in the case of metastable states. A celebrated 
example where this description has proved extremely useful is the two-states kaon com- 
plex. If this system is influenced by an external interaction, the Hamiltonian operator 
can be written as 

H = H +V (19) 

where TCo and the perturbation V are, in general, two non Hermitian and noncommuting 
operators. Thus, we can apply the formulas of the previous section to expand U(t) in 
terms of the eigenf unctions |Ks), |Kl) of Ho 

Ho|Ks)=As|K s ) 

(20) 

n \K L ) =A L |K L ) 
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where |Ks) and |Kl) are the right eigenvectors. The same matrix operator Ho has also 
two left eigenvectors with the same eigenvalues 

(K' S |W =A S (K^| 
(K^ =A L (KiJ . 
The set {|K' S ), |Kl)} is the reciprocal set of {|Ks), |Kl)} in the sense that 



then 



|K' S ) = J_ (|Ks>-x|K L » 

Vi-lxl 2 

\K)= (|K L )-xlKs)) 

Vi-M 2 



where 



and 
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' (|K' S )(K S | + |KL)(K L |) 



i-lxl 2 



(21) 



(K' S |K L ) = = (K' L |K S ) . (22) 
If we normalize all the eigenvectors to 1 and denote the overlap as 

X = (K L |K S ) (23) 



(24) 



(K / S |K S > = ^1-| X | 2 = (KL|K L > (25) 



(K' L |K' S ) = - X . (26) 
Therefore, we have the following equivalent decomposition of unity 

J=|K s )(K s | + |Ki j )(Ki j | = 
= |K L )(K L | + |K' S >(K' S | = 



1 (27) 
(|K S )(K' S | + |K L )(K^|) = 



In the particular case of the evolution operator, the generalized Cauchy's formula is 
written in the form 



U{t) = ^ j e- izt Q{z) dz 



(28) 



where the resolvent operator is 



l 



l 



z X - H zl- [Ho + V] 



(29) 



and T is a closed curve encircling all the complex eigenvalues of the total Hamiltonian 

n. if 

i 



(30) 



zl-Ho 

is the resolvent operator of Hq, an analogous expansion to Eq.(6) for Q(z) holds by 
substituting Ho and V for the operators A and B respectively. Now, the matrix elements 
will be obtained by 

U a p = (K' a \exp[-iHt]\Kp) (31) 
where the greek letters a and (3 are fixed and a, (3 e {S, L}. From Eq.(12) we have 



(K^|exp[-iHt]|K > = 



E E «\ V V fe^M - fe^M v|K,) x 



n=0 { M ,n} 



(K^|K M1 ) (K^ 2 |K M2 ) - (K^JK^; 



x 



2^ / dze-** [(z - K)-\z - \^)-\..{z - X^y'iz - Xp)~ x ] = 



(32) 



V|K^)x 



n=0{ M ,n} v Mil ^1/ \ A* 



X 



n-1 



7e(A Q ) + ^7e(A M J + ^(A^ 



I/=l 



Here {/i, n} = {/zi, /Z2, A*n-i} an d the indices are varying in the set {S, L}. We 
can give a more explicit expression to the quantity (K^,| exp[— iHt] jK^) with the use of 



formulas (23), (24), (25) 

U a p = {K , a \exp[-iHt]\Kp) = 



E 

n=0 



i-lxl' 



(n-1) 



(K'JVIIKsXK'sI + IKlXK'JJV 



(33) 



[ IKsXK'sI + |K L >(K L | ]...[ |Ks>(K' s | + |K L )(K^| ] V^) x 



n-l 



X 



Being A a , A^ and A^ [y e {1, 2, n — 1}) equal to As or Al, Eq. (33)will be therefore 
rewritten as 



i 



n=0 



yji-\xf 



TT ^ Ty (K^|V[|K s )(K' s | + |K L )(Ki j |]V 



[ |Ks)(K s | + |K L )(K' L | ]...[ |Ks)(K' s | + |K L >(K L | }V \K P )[K(\ S ) + K(\ L ) ] 



where TZ(Xi) is the residue at z = Aj, Aj G {As, Al}, of the function 



F{z)=e- lzt [{z-\s) r {z-\^) s ] 



-l 



(34) 



(35) 



and r and s are positive integer numbers subject to the condition r + s = n + 1. In a 
general theory 



trH T ^[trHo] 2 -MetHo 
a s,l = r (36) 



so As 7^ Al and all eigenvalues are not degenerate. In this way if n > 2, As and Al are 
not simple poles for F(z) and a direct calculation of residues gives 

1 cT" 1 



and 



K(\ L ) = 



(r - 1)! dz r ~ 

1 d 5 " 1 
(s - 1)! dz 6 



T [e-^(z-X s )- s } 



T [e—\z-\^] 



z=\ s 



z=\t_ 



(37) 



(38) 
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We display the use of these formulas to calculate the second order term in expansion in 
the particular case a = S and j3 = L 



U, 



(2) 



S,L 



(K^| V |K s >(K's| V |K L ) [K {2) (\ s ) + ft ( i)(A L ) ] + 

'i-lxf 

+(K^| V |K L )(K^| V |K L ) [ft (1) (A s ) + n (2) (X L ) [ 



(K' S |V|K S )(K^|V|K L ) 



i-ixr 

+ e-^(A L -A s )- 2 



-iA s t 



[(Al-As)" 1 + (Al-As)" 2 ] + 



+ 



1 



(K^|V|K L )(K^|V|K L ) 



i-lxl 

+ e- AB '(A s -A L )- 



-iA L t 



[(As-Al)" 1 + (A S -A L )- 2 ] + 



(39) 



Here we have labelled 7£(Ai) with the subscripts (1) and (2) to stress that Aj is a 
first or a second order pole. Obviously, the presence of higher order terms involves 
higher order poles, but, in the case of the neutral K meson system, the expression of 
U a [3 is not so cumbersome as in the general case. In view of this consideration, this 
perturbative approach is extremely useful in the description of evolution of the two- 
states kaon complex [4]. But, as mentioned in [ 1 ], the use of the complex analysis is 
more convenient and it can be successfully applied also when the usual time-dependent 
perturbation theory fails. 



III. Concluding Remarks 



Outside the realm of particle physics, there are many other cases of unstable systems 
influenced by external interactions, where the previous approach becomes indispensible. 
For example, in modern quantum optics, it seems particularly important to analyze 
the (para) magnetic resonance [5], and in general to describe the two states (spin-up, 
spin-down) involving electrons and protons with dissipation. Presently, it provides 
the theoretical framework to study a multitude of effects involving laser dynamics. 
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Nevertheless, unstable two level systems in interaction with other degrees of freedom 
require the strategy outlined before. The system is, in fact, an open system and its 
dynamical behaviour under the influence of an external interaction can be described 
only redefining the evolution operator 

OitYAifiO- 1 ^) = U{t) = exp[-itH/h] = e A+B . (40) 

The transformation 0(t) decouples the internal degrees of freedom from the motion 
of the center of mass and provides a time independent Hamiltonian H = OHO -1 — 
%hOO~ x . The operators A and B are introduced to clarify the mathematical structure 
of the calculation below. The new time evolution operator li{t) is then determined 
by the exponential factorization of two non commuting (sometimes non Hermitian) 
operators. It is evident now, the importance of the method outlined above which turns 
out to be particularly compelling as far as the physical interpretation is concerned. It is 
worth discussing in connection with the algebraic approach. This method makes use of 
the parametric differentiation of the exponential of an operator and of the commutation 
relations in the context of the Baker-Campbell-Hausdorff (BCH) formula [6]. Then we 
can think to separate the center of mass A part evolution to factorize U(i) according to 

U(t) = e A W(t) . (41) 

Thus the complete time evolution of the two-level system is then based on the remain- 
ing determination of the operator W(t) which contains the influence of the external 
interactions on the internal dynamics. On this point, to work out W(t) of Eq. (41) we 
consider the operator 

Q(X) = exp[A(A + B)] = e XA W(X) , (42) 

and restrict to A = 1 at the end. Differentiation of Eq. (42) with respect to A leads to 
the differential equation 

^ = (e~ XA B e XA ) W(X) ~ (B - X[A, B]) W(X) (43) 

with the initial condition W(X = 0) = 1 and using the identity 

oo 

e -XA B e XA = >^K n , K = B , K n+1 = [K n , A] (44) 

n=0 
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which holds for any two operators A, B. Under general assumptions, it may be written 
as a matrix equation of the form, 

fdWn/dX dW 12 /dX\ ( P-X[Q,P] R \ / W u W 12 \ , . 

\dW 21 /dX dW 22 /d\) \ R -P + X[Q,P]J \W 21 W 22 ) ■ { 0) 

This result assumes that the factorization is able to select the non vanishing commutator 
[A, B] = [Q,P], that is, however, a c-number, whereas R is the remaining part in B. 
An exemplification is represented by 

A =Q1 

(46) 

B =Pa 3 + Ra t . 

Eq. (45) is an operator-valued system of differential equations, but it contains only 
commuting operators so that we can treat it as an ordinary differential equation with 
the initial condition Wij(X = 0) = 1. 

Inserting the equation for dWn/dX into the equation for dW 2 \/dX and similarly 
for dW 22 /dX into that for dWi 2 /dX, one gets 

%^={ J R 2 -[Q,P] + (A[Q,P]-P) 2 }^ii , N 

(47) 

%^ ={R 2 + [Q, P] + (A [Q, P] - Pf}W 22 . 
After the introduction of the parameter 

9 = 2[ftp] ( 48 ) 
and the change of the variable y = (X[Q, P] — P)^/2/[Q, P], Eq. (47) becomes 

(49) 

with the initial conditions Wn(X = 0) = W 22 (X = 0) = 1 and 

P 



dWn 
dy 



dW 22 



dy 

A=0 



IQ,P] ■ 

A=0 



(50) 



The solution of Eqs. (49) is a linear combination of parabolic cylinder functions. The 
total operator U{t) is given by 

«W-pw(S wll) ■ ( 51 ) 
li 



It remains to cancel the initial unitary transformation 0(t) in Eq. (40) to obtain the 
exact expression for the time evolution operator: 



Uit) = O' 1 (t) U{t) 0(0) = exp {A} O' 1 



(t)W(t)O(0). 



(52) 



An instructive consistency check is to turn off the external interaction by setting 9 
to zero. In this case the stable two-level system should be recovered. 6 = implies 
immediately 



as it was to be expected. The first term describes the free motion of the system, and the 
second term contains the internal transitions. It is also possible to derive an expansion 
for small 9. But the result is then difficult to understand since it contains various 
combinations of error functions. 

In this paper, we have analyzed the dynamical evolution of unstable systems under the 
influence of external interactions. The generalization of the complex spectral theory is 
proposed to account for these unstable open systems. Furthermore, the results of the 
coupled dynamics of the internal transitions and the center of mass motion are worked 
out with the algebraic expansion of the time evolution operator. 



One of us (M. V.) wishes to thank the warm hospitality of the Department of Mathe- 
matics of the University of Basilicata, Potenza. 



U(t) = exp {— itHcm} exp {— itV} 



(53) 
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